Order-parameter fluctuations in the frustrated Heisenberg model on the square lattice 
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The T — dynamics of the two-dimensional s — 1/2 Heisenberg model with competing nearest- 
neighbor (Ji) and next-nearest-neighbor (J2) interactions is explored via the recursion method, 
specifically the frequency-dependent fluctuations of the order parameters associated with some of the 
known or suspected ordering tendencies in this system, i.e. Neel, collinear, dimer, and chiral order. 
The results for the dynamic structure factors of the respective fluctuation operators show a strong 
indication of collinear order at J2/J1 > 0.6 and a potential for dimer order at 0.5 < J2/J1 ^ 0.6, 
whereas the chiral ordering tendency is observed to be considerably weaker. 
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The impact of quantum fluctuations on the zero- 
temperature phase diagram of a quantum many-body 
system tends to be strongest if that system contains 
competing interactions. A particular, sometimes exotic 
phase may be stabilized by quantum fluctuations in the 
presence of interactions that frustrate each other's order- 
ing tendencies. The two-dimensional (2D) Heisenberg 
antiferromagnet with nearest-neighbor (mi) and next- 
nearest-neighbor (nnn) coupling on the square lattice, 
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has been a prominent object of study in this contextE 

The Neel long-range order (LRO) present in the ground 
state of the nn model disappears at some critical cou- 
pling ratio, J2/J1 — 0.4, and is replaced, at sufficiently 
large J2/J1, by Neel LRO within each of the two nnn 
sublattices. The latter is preceeded, at J2/J1 ^ 0.65, 
by collincar LRO, which breaks the (discrete) rotational 
symmetry of H on the lattice but not yet its (continu- 
ous) rotational symmetry in spin space. The Neel and 
collinear order parameters (OPs) are described by the 
operators, 
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At intermediate coupling ratios, 0.4 < J2/J1 0.65, 
the Neel and collinear ordering tendencies keep each 
other at bay and thus make the frustrated ground state 
susceptible to different kinds of ordering potentialities. 
Dimer order p u twist order ,□ and chiral orderH have been 
proposed in this context. The dimer and chiral OPs con- 
sidered here are defined by the following expressions in 
terms of spin operators: 
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where D T — S r -S r _|_ Xj Xr — Zr,r+x,r+x+y — ■Zr.r+x+y.r+y , 

^i.j, k = *(SfSr - Sr S + + S+S w - Sr S + + S+Sr - 
S^S^~). They probe the long-range phase coherence 
of singlets stacked in columns along the y-axis and 
the handedness of the spin configuration on a plaque- 
tte, respectively, in the ground-state wave function \G). 
Whether any one of these types of LRO is, in fact, real- 
ized, or whether a ground state with shorL-xange correla- 
tions of the resonating- valence-bond type,Bl3 for example, 
is stabilized, has not been determined for certain. 

The absence or presence of a specific type of LRO de- 
termines whether the associated OP correlation function 
decays to zero or not. In the finite-cluster data, the rele- 
vant information on the asymptotic behavior is best cap- 
tured by the expectation value of the squared OP. The 



problem here is to find a meaningful reference point for 
any enhancement in that quantity.!] 

An alternative avenue to comparing the different or- 
dering tendencies in the spin-frustrated ground state of 
([!]) is to explore the dynamic (i.e. frequency-dependenid 
fluctuations of any proposed OP. The recursion methodS 
in conjunction with recerdjhz. developed techniques of 
continued-fraction analysisoEJ is very suitable for that 
purpose. Here the dynamical information is derived from 
the finite-size ground-state wave function. No excited 
states have to be computed. This is an important advan- 
tage for the study of systems with complicated spectra 
and with potential OPs that have widely varying sym- 
metry properties, as is the case here. 

We investigate the fluctuations of the four OPs defined 
in (^]) and (|^) as they manifest themselves in the dynamic 
structure factors 
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where F^ = iV _1 / 2 E r e* q ' r A r is the fluctuation operator 
associated with a given OP, and A T stands for N T , C r , 
D r , or Xr, as defined above. 

The recursion algorithm in the present context is 
based on an orthogonal expansion of the wave function 
l^q W) = It produces (after some interme- 

diate steps) a sequence of continued-fraction coefficients 
A^(q), A^(q), . . . for the relaxation function, 
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which is the Laplace transform of the symmetrized cor- 
relation function K(i^(t)F^ q )/(F^Ff q ). The T = 
dynamic structure factor is then obtained from (|^) via 

S AA (q, w) = 4(^4)6(0,) lim &[c£ A (q, e - ioj)} , (6) 

where(F q 4 F^ 4 q ) = Saa(q) is the integrated intensity.0 

All the results presented here are for clusters of N = 
4x4 sites with periodic boundary conditions. Extreme 
care must be exercised in separating finite-size effects 
from properties that reflect the physics of the infinite sys- 
tem. This distinction can be made with more confidence 
for coupling ratios J2/J1 ^ 0.7 than at J2/J1 > 0.7, 
where the gradual decoupling of the two nnn sublattices 
causes a crossover in the finite-size effects. 

The Neel OP fluctuations are probed by the dynamic 
structure factor SjVAr(q, oj) at the wave vector q = (ir, 7r). 
This quantity is shown in Fig. 1 for various coupling ra- 
tios. The presence of Neel LRO at J2/J1 = implies 
that S n n (ft , ft , is governed by a zero-frequency peak. 
Quantum fluctuations split the ground-state level for fi- 
nite N. In a 4 x 4 cluster the finite-size gap is known 
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to be of magnitude AE/Ji ~ 0.57.t2l The peak position 
of the curve for J2I J\ = must be interpreted with this 
fact in mind. 

What are the effects of the rain coupling on the Neel 
OP fluctuations? For < J2/J1 i5 0.4, i.e. over the es- 
timated range of the Neel phase, we observe only small 
changes in peak position and line shape. Then the Neel 
OP fluctuations begin to change rapidly in two stages: 
(i) Over the range 0.4 < J2/J1 ^ 0.6, the peak posi- 
tion moves to higher frequencies at an accelerated rate, 
the linewidth shrinks, and the integrated intensity (not 
shown) drops to 32% of its value at J2/J1 = 0. This 
signals the presence of some non-Neel type ordering ten- 
dency which supports well-defined Neel modes at increas- 
ingly high frequencies, (ii) At J2/J1 ^ 0.6 the linewidth 
of Snn{^, tt, k-0 grows rapidly, while the peak position 
moves further up and the integrated intensity continues 
to fade away quickly. As the two nnn sublattices begin to 
decouple, the system ceases to support well-defined Neel 
modes. 

The dynamic structure factor <S'cc(0, 0, w), which de- 
scribes the collinear OP fluctuations, is shown in Fig. 2. 
At J2/J1 = 0, we observe a fairly sharp collinear mode 
at uj/Ji ~ 3.0. As the Neel ordering tendency weakens 
with increasing J2/J1, the collinear mode shifts to lower 
frequencies, while its line shape broadens considerably. 
The width reaches a maximum at J2/J1 — 0.4. Between 
here and J2/J1 — 0.55, where the competing dimer and 
chiral ordering tendencies are at their peak, the collinear 
mode moves to u> = 0, and the integrated intensity more 
than triples in relation to its value at J2/ J\ = 0. In the 
interval 0.55 < J 2 /J\ < 0.7, the function S C c(0,0, a;) 
transforms into a narrow central peak, and Scc^O) in- 
creases by another factor of ~ 1.6. This clearly reflects 
the onset of collinear LRO in the infinite system. 

The dimer OP fluctuations as described by the dy- 
namic structure factor Sdd{^, 0,o;) and shown in Fig. 3 
resemble those of the collinear OP with respect to line 
shape and peak position for as long as the Neel ordering 
tendency is perceptible in the ground state ( J2/J1 ^ 0.4). 
Both modes become soft and very broad at J2/J1 — 0.55, 
but then they part company. While the collinear mode 
has been observed to transform into a high-intensity nar- 
row central peak, the dimer mode, which has reached 
its maximum intensity here (~ 1.8 times its value at 
J2I J\ — 0) , broadens further and and loses intensity very 
rapidly. It literally dissolves as the inter-sublattice cor- 
relations begin to weaken at J2/J1 ^ 0.7. Nevertheless, 
the softness of the dimer OP fluctuations at J2/J1 — 0.55 
in the 4x4 cluster is consistent with dimer LRO in the 
infinite system. 

The chiral OP fluctuations are based more significantly 
on intra-sublattice correlations than the dimer OP fluc- 
tuations and, therefore, evolve differently. This is illus- 
trated in Fig. 4. At J2/J1 = the dynamic structure 
factor S xx (0,0,uj) exhibits a sharp mode at lo/J\ ~ 3.8. 
At J2/J1 — 0.55, the peak position has moved down 
to lu/Ji ~ 1.6, while the linewidth has increased only 



slightly, and the integrated intensity has grown to a max- 
imum value of ~ 2.3 times its value at J2/J1 = 0. Then 
the peak position starts to move back out to higher fre- 
quencies, the line shape begins to broaden, but less so 
compared to that of the dimer fluctuations, and the in- 
tegrated intensity drops rapidly. The minimum gap of 
the chiral mode is perhaps too large to be entirely at- 
tributable to a finite-size effect, which would indicate that 
the observed chiral ordering tendency does not turn into 
chiral LRO as N — > 00. 
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FIG. 1. T = dynamic structure factor (|4|) normalized 
by its integrated intensity for the Neel fluctuation operator 
Fq at q = (tt,ty) of the Hamiltonian ([l]) at various values of 
the coupling ratio J2/J1, obtained via strong-coupling con- 
tinued-fraction reconstruction from the coefficients Ai, As 
and a Gaussian terminator as explained in Refs. pj,[L0[ 

FIG. 2. T = dynamic structure factor (^) normalized by 
its integrated intensity for the collinear fluctuation operator 
at q = (0, 0) of the Hamiltonian (|l|) at various values of 
the coupling ratio J2/J1, obtained via strong-coupling con- 
tinued-fraction reconstruction from the coefficients Ai, As 
and a Gaussian terminator. 



FIG. 3. T — dynamic structure factor (|J) normalized 
by its integrated intensity for the dimcr fluctuation operator 
F^ at q = (it, 0) of the Hamiltonian (|l|) at various values of 
the coupling ratio J2/J1, obtained via strong-coupling con- 
tinued-fraction reconstruction from the coefficients Ai, As 
and a Gaussian terminator. 



FIG. 4. T — dynamic structure factor (|4|) normalized 
by its integrated intensity for the chiral fluctuation operator 
Aq at q = (0, 0) of the Hamiltonian (^) at various values of 
the coupling ratio J2/J1, obtained via strong-coupling con- 
tinued-fraction reconstruction from the coefficients Ai, As 
and a Gaussian terminator. 
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